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Abstract 

Pupil imaging performance is analyzed from the perspective of physical optics. A multi-plane diffraction model is 
constructed by propagating the scalar electromagnetic field, surface by surface, along the optical path comprising the 
pupil imaging optical system. Modeling results are compared with pupil images collected in the laboratory. The 
experimental setup, although generic for pupil imaging systems in general, has application to the James Webb Space 
Telescope (JWST) optical system characterization where the pupil images are used as a constraint to the wavefront 
sensing and control process. Practical design considerations follow from the diffraction modeling which are discussed in 
the context of the JWST Observatory. 
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1. INTRODUCTION 

A pupil imaging optical sub-system is designed to image the aperture stop of the parent optical system. The exit pupil 
(hereafter referred to as the “pupil”) is defined as the image of the aperture stop as seen from the detector. In practice, 
the pupil image is an important diagnostic and systems engineering tool for assessing optical system performance. 
Factors affecting optical performance which are readily obtained from pupil image include: pupil geometry, obscuration 
profiles, radiometric data, and optical sub-system alignment. Obtaining this information from the pupil image will be 
limited by the pupil image quality which can be quantified and toleranced using conventional ray-trace techniques. But 
just as important in this assessment of image quality is predicting how image detail produced by the pupil imaging 
system will be affected by image ringing and other diffraction effects, recalling that diffraction occurs at the system stop, 
the object being imaged. Therefore, designing a pupil imaging system should not just include the standard ray-trace 
analysis, but also a multi-plane diffraction model. Modeling these effects falls under the domain of physical optics. 

In this paper 1 we discuss pupil imaging performance from the perspective of physical optics. In Section II, a wave- 
optics model is constructed by propagating the scalar electromagnetic field along the optical path comprising the pupil 
imaging system. In Section III, the model results are compared with pupil images collected in the laboratory. The 
experimental setup, although generic for pupil imaging systems in general, has application to the James Webb Space 
Telescope (JWST) 2 optical system where the pupil images are used as a constraint to the wavefront sensing and control 
process. 3 


2. WAVE-OPTICS MODEL 

The optical system used for the analysis is illustrated in Fig. 1. The aperture stop in this system is the primary mirror. A 
wave-optics model is constructed by propagating the scalar electromagnetic field from surface to surface. For clarity, 
the PIL (pupil imaging lens) is excluded from Fig. 1 to emphasize that the complex amplitude of the scalar field at the 
focal plane of the mirror is given by the Fraunhofer propagation of the complex pupil, A(x,y)e lk<p(x,y) : 

u = 3{A(x,y)e ik<p(x ' y) }, (1) 

where (x,y) denote pupil coordinates and k = 2nl X for monochromatic imaging. The complex pupil function consists 
of the binary aperture function, A(x,y ) , and a phase term, <p(x,y ) , representing the optical path difference for the case 
of an aberrated mirror surface. 
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Fig. 1. Optical System 

To find the field after propagation through the PIL, we begin first by specifying the field, u x , immediately incident upon 

the first PIL surface 1 as illustrated in Fig. 2. For simplicity, the PIL is shown as a singlet in Fig. 2 but in 
implementation the PIL was actually a doublet from Edmund Scientific. 4 



Fig. 2. Optical System with the Pupil Imaging Lens (PIL) 


As shown in Fig. 2, the field immediately preceding the first PIL surface is at a defocused plane, or equivalently, in the 
near field or Fresnel zone (this correspondence is discussed in detail elsewhere 5 ). The field, u x , at the defocused plane- 
1 is 


where, 


u x =Z{A(x,y)e ika <*■»}, 
9{x,y) = (p + cp d , 


( 2 ) 

( 3 ) 


and <p d (x,y) = a d (x 2 + y 2 ) is the quadratic (defocus) phase term with a d in units of X . We note that it is also possible 
to place m, much farther into the near field, i.e., closer to the aperture stop, by using the Huygens-Fresnel propagator (a 
derivation in the pupil imaging context is given elsewhere; 6 ® denotes convolution): 

«, = 3{A(x,y)e lklf(x - y) } ® e^ +/) . (4) 
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However, expression (2) is the most relevant for our model since the experimental setup has the PIL placed adjacent 
from best focus. We make the point that after inclusion of the PIL lens in Fig. 2, the PIL is under-filled, and therefore, 
the primary mirror remains the aperture stop of the optical system. 


To proceed, the total phase delay imparted onto the wavefront upon propagation through the PIL, from Fig. 3 is 

© = nt(x,y) +t 0 - t(x,y ) + y(x,y) 


= t 0 +(n - 1 )t{x,y) + y(x,y), 

where nt(x,y ) is the phase delay due to the lens with index of refraction, «, and 
t 0 -t(x,y) is the phase delay due to the intervening space from both the front and back 
of the lens. The y(x,y) term represents the overall contribution due to the optical path 
difference as a result of the PIL aberrations. Therefore, the field at surface 2 of the PIL 
shown in Fig. 2 is given by the product: 

u 2 =u x e ike{x ' y) = u x (e iktQ )(e ik{n ~ mx ' y) ) e ikr{x ' y) . (6) 

The first exponential in (6) is a constant; the second term is the phase delay due to the 
lens with index of refraction, «, and the third term is the phase delay due to the lens 
aberrations. 



( 5 ) 


To continue, we must express t(x,y) in terms of the pupil coordinates, (*,>>)• A 
standard derivation 7 follows from the construction of Fig. 3. For completeness, we 
include a derivation of t(x,y) in the Appendix. Using these results, Equation (B30), 

t(x,y) is expressed: 


t(x,y)=t 0 + 


(x 2 +y 2 ) 


_1 1_ 

R 2 R , 




The second exponential term of (6) is thus 


ik(n - 1 )t(x, y) = ik(n- 1)/ 0 + ik(n- 1) 


1 1 


\* 2 +/) 


^2 ^1 J 


but note that the focal length of the lens is given by 

1 


/ 


= -("-!) 


f \ 1 ^ 

\^2 ^1 J 


which allows (8) to be expressed in the simpler form: 


ik(n-l)t(x,y) = iknt 0 - ikt 0 -ik 


(* 2 +/) 

2/ 


Therefore, the field at surface 2 of the PIL as shown in Fig. 2 is given by the product: 


-/ — (or +> z ) 

o ikt »\(o- ikt »\(o- iknt »\(0 2 f ) e iky(x,y) 


u 2 = u \ ( eikt ° )( e ~ ikt ° )(e iknt ° )(e 


Fig. 3. Phase Delay due 
to the PIL. 


( 7 ) 


-i-L( X 2 + y 2 ) 

= 21 e iky(x ' y) . 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 


In the back focal plane of the lens the field, w 3 , is obtained by propagation of u 2 using the Huygens-Fresnel integral: 


w 3 = u 2 (8) 


i-r-(x z +y l ) 

o 2z 


= —e ib j jdx'dy’[u 2 ]e 


i—[( X - X y + (y-y’) 2 ] 


= constant x J J dx'dy'[u 2 \ e 


'Yj( x ' 2+y ' 2 ) -i^(xx'+yy') 


( 12 ) 


e 2 


where the second form is obtained after expanding the exponential in the preceding integral. Substituting u 2 from (11) 
and noting that in the back focal plane, z -» / , the field at w 3 is then 
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( 13 ) 




w 3 « J | dx'dy’\u x ( e' knt ° ) e 2/ 

-00 

00 

= e lk "'° | \dx'dy'[u x e ikrtx ’ y) ]e 


.. . . 'Trf *' 2 + y' 2 ) -j— (xx'+yy‘) 

e -krU.y)- i g 2/ g / 
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where the quadratic phase factors cancel and the constant term, e * 
this result, the field at w 3 is then expressed: 


, is factored outside of the integral. Normalizing 


u i = \\dx'dy'[u x e lkr(x ' y) ]e 


-/-(xx'+^y') 


(14) 


which is a Fourier transform: 

m 3 =3{m 1 ^ m } = 3{« 1 }03^ (15) 

expressed on the RHS as a convolution product of the Fourier transform of the incident field and the Fourier transform 
of the PIL complex amplitude defined by the lens aberrations. Substituting from (2) for u x we get the result 

M 3 =3{3[^(x,^)e iW( ^ ) ]}<8)3{e i * /(x>>) }. (16) 

The result in (16) can be further simplified by noting the mathematical identity associated with the Fourier transform: 

g(-v) = 3{3[g(x)]}, (17) 

i.e., the Fourier transform of the Fourier transform gives back a scaled and rotated version of the original function. 
Applying (17) to (16) gives the pupil image complex amplitude at the detector, plane 3: 


/ 


= u } =[A(-v x - v y )e‘ kS{ - v ”- v ’ ) ]®Z{e ikr(x - y) } . 


(18) 


1 pupil ,ca ** 3 L-^v w y > 

The result in (18) states that the pupil image intensity is given by the product of the parent optical system pupil complex 
amplitude, pupil ca optical system , convolved with the PIL amplitude spread function, ASFp IL = 3{e lkr(x ' y) } , i.e., the PIL 
amplitude spread function without a PIL aperture response: 

1 pupil, , mg = 1 pupil, cu ■ Tpupn.cc = II PW‘L ® ASF p,l II 2 • ( 1 9) 

The prime on ASFp IL , denotes that it is not the true complex amplitude of the PIL since no aperture term for the PIL is 
included in the right-hand term of (18). The reason is that the PIL imaging was assumed to be under-filled as it is for the 
JWST PIL optical design, and the laboratory system we have constructed in Fig. 2. However, there is no more difficulty 
in allowing the PIL to be overfilled, in this case Equation (18) then becomes: 

I pupil, ca =«3 = [A(-v x , -V y )e ik ^ •-^ ) ] 0 Z{A p[L (x, y)e ikr(x ' y> } . (20) 

An important limit follows from (18) for the case that the PIL lens has no aberrations: the pupil image complex 
amplitude is simply a scaled version of the complex aperture stop: 

Ipu Pi ,.cc=A-v x -v y )e ikS ^- vA . (21) 

Therefore, the pupil image collected by the detector at the focal plane follows from the complex conjugate square of (18) 


I pupil, ca I pupil, ca 


pupil, img 


= A(-v x ,-v y f 


( 22 ) 


showing that when no aberrations are introduced by the PIL (y(x,y) = 0 ), the image is essentially the image predicted 
by geometric optics, i.e., the pupil image is a scaled version of the aperture stop. 


This conclusion is essentially a well-known result 7 (p. 130, Equation (6-3)) and forms the basis of Fourier optics, since 
this result allows a system to be compressed from a complex multi-plane system to a simple system consisting of 
entrance and exit pupils and focal planes, i.e., the optical system may be regarded as a black-box. Note especially that 
for the case of a perfect PIL lens (y(x,y) = 0 ), the primary mirror aberrations, <p(x,y) (see Equation (3)), do not affect 
the pupil image quality. This is due to a cancellation of the complex phase term in the (19). The result is perfectly 
reasonable to within the paraxial assumption since in essence, the PIL is imaging the primary mirror boundary which is 
the aperture stop in the system. 
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3. EXAMPLES 

Given the result in (18), it straightforward to calculate the image quality for an arbitrary PIL system given the optical 
path difference induced by the PIL lens itself, y(x,y), where y(x,y) is available from a ray-trace model or phase- 
retrieval measurement. 3 For the system modeled in Fig. 2, our ray-trace model predicts y(x,y)& 1 wave of spherical 
aberration where lwave = 670 nm. A comparison of both the model and data are shown in Fig. 4 where good qualitative 
agreement is shown upon comparison of the data and model. No attempt was made to account for detector and sampling 
effects in the correlation of data and model. 



Fig. 4. Comparison of Data (a) and Model (b). 


It is important to point out that our coherent imaging model predicts the well-known ringing 8 in the shadow, although 
detector sampling in our laboratory measurements did not permit it measurement. This is illustrated across one of the 
segment edges in the JWST model in Fig. 5. 



Fig. 5. Coherent ringing across a segment edge in the JWST model. 








As a final application of the wave-optics pupil imaging analysis, we calculate the predicted image-quality of the JWST 
PIL design. 9 To determine the in-situ PIL optical path difference, y(x ,y ) , we use a full-ray-trace model of the JWST 
flight design 10 including the Near-Infrared Camera (NIRCam) 11 in combination with the PIL. The y(x,y) is taken to 
be the difference between the exit pupil optical path difference (OPD), first with the PIL in the design, and then with the 
PIL removed from the design. The results are shown in Fig. 6 in the Malacara basis set. 12 


y = (opd in - opd out) Wavefront from Fit: 
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Fig. 6. JWST PIL optical path difference, y(x,y) . 


The numerical values of the decomposition shown in Fig. 6 are listed in Table 1 in a variety of units. 
Table 1. Zemike decomposition of the JWST PIL y(x,y) 


Z4 

Astigmatism 
(1st order 45°) 

Z6 

Astigmatism 
(1st order 0°) 

Z7 

Trefoil 30° 

Z8 

Coma Y 

Z9 

Coma X 

Z10 

Trefoil 0° 

UNITS 

-2.1715 

3.9015 

0.0594 

-0.4861 

-0.3603 

0.0961 

pim 

-1.1524 

2.0705 

0.0315 

-0.2580 

-0.1912 

0.0510 

waves (2 = 
1.884 fxm) 

-2171 

3901 

59.42 

-486.15 

-360.27 

96.06 

nanometers 


As shown in Table 1 and Fig. 6, the dominant contributor to y(x,y) is astigmatism. Substituting y(x,y) into (18) and 
then using the 18-segment mask of the JWST design for A(-v x ,-v y ) , a “cross-stitch” ringing pattern is induced from 
the astigmatic nature of y(x,y ) . The result is shown in Fig. 9. The structure of this astigmatic ringing pattern is better 
understood by using a circular aperture mask for the aperture stop, but with the same y(x,y ) . This result is shown in 
Fig. 8, where it is apparent that the astigmatism introduces the rather unusual “cross-stitch” ringing pattern. 

4. SUMMARY AND DISCUSSION 

We have presented a paraxial wave-optics model of coherent pupil imaging and have compared modeling results with 
pupil images collected in the laboratory. An interesting result is that coherent ringing may be traced to the PIL lens 
aberrations, or optical path delay. Indeed, for the limiting case that y(x,y) = 0 , the resulting image produced by the PIL 
is essentially the image predicted by geometric optics, i.e., the pupil image is a scaled version of the aperture stop. 

The practical implications are that the pupil image will not be significantly degraded by diffraction effects if two 
conditions are met: optical aberrations are well-corrected and the subsequent optical system apertures do not clip the 
wavefront coming from the mirror’s surface, recalling the assumption that the PIL lens itself was assumed to under- 
filled. 13 Therefore, the pupil image that is produced from a well-corrected PIL with sufficiently oversized apertures will 
closely approximate the pupil image predicted by geometrical optics, for obscurations located on or near to the system 
aperture stop. 
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The experimental setup described in this paper, although generic for pupil imaging systems in general, also has 
application to the James Webb Space Telescope (JWST) optical system characterization where the pupil images are used 
as a constraint to the wavefront sensing and control process. We have determined the in-situ JWST PIL optical path 
difference using the flight ray-trace model and have used these results to predict the exact nature of the coherent ringing 
produced by the PIL imaging system. An obvious concern due to the ringing in the PIL images is the ability to resolve 
or register pupil edges. Although is not necessarily a concern for systems with a simple circular aperture stop, this is 
definitely a concern for complex apertures such as the JWST 18 segment aperture. We note however that using a 
straightforward image processing technique such as deconvolution can reduce these ringing effects. An example of this 
post-processing step is shown in Fig. 7. 


pypii edge data: deconvolution results: 



Fig. 7. Pupil image coherent ringing reduction using one deconvolution iteration. 

A much simpler and more effective technique for reducing the edge ringing is to image the pupil in broadband light. In 
this case the temporal coherence length, L c : 14 

L c =A 2 / AA, (23) 

tends to zero and the coherent ringing disappears from the images. Given that the signal-noise-ratio of the broadband 
light will be greater than in the narrow-band case, the resulting pupil imaging exposure times can be significantly 
reduced as well, at the expense of additional chromatic imaging effects. 
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Fig. 8. Image Predicted by Coherent Model of JWST Pupil Imaging. 
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APPENDIX 

For completeness, a derivation of t(x,y) is given. Other derivations may be found in the literature. 7 To begin, we adopt 
the convention shown in Fig. 10(a) that the radius of curvature is positive for the first surface (convex) and that R< 0 
for the second surface (concave). To find the explicit form of t(x,y ) in terms of the pupil coordinates, ( x , y ), the 
geometry of Fig. 3 is further divided into three regions as shown in Fig. 10(b). 




(a) 



Fig. 10. (a) Radius of Curvature Sign Convention; (b) Regions of Phase Delay due to the PIL. 
Therefore, we consider 

t(x,y) = t x +t 2 +t 3 (B24) 

and 


~ ^oi ^02 ^03 j (B25) 

where t 0 defines the thickness of the lens long the optical axis. Fig. 11 shows additional detail on Region 1 of Fig. 
10(b). 


8 





Fig. 11. Detail on Region 1. 


From Fig. 11 the explicit form of t x (x,y) is obtained: 

= oi —[^i (-^5 jv )] 

= *0, - [* , - n /^ 1 -(x 2 +/)] (B26) 

=? 01 -/? 1 [l- A /l-(x 2 +/)/*?]. 

The expression for t x (x,y) is further simplified using the binomial approximation which is equivalent to the small angle 
approximation of paraxial optics: 

hix,y) = *oi - *, [l ■ - (1 ■ - (x 2 + y 2 ) / 2R] )] 

= t m -R,[{x 1 +y 1 )l2R])\ (B27) 

= t m -(x 2 +y 2 )/ 2R X . 

The expression for t 2 (x,y ) is similarly obtained: 

t 2 (x,y)=t 02 +(x 2 +y 2 )/2R 2 , (B28) 

after factoring out a (-\)R 2 term consistent with the sign convention of Fig. 10. Substituting (B27) and (B28) into 
(B24) and then using the fact that t 2 = t 02 gives 


t(x, yj — t 0i + t 02 + 1 0 3 + 


{x 1 +y 2 )[ 1 1 


R 2 R , 


and after using (B25) the result simplifies to 


(x 2 +y 2 ) 1 1 


2 I R 2 R t 


(B30) 
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